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ABSTRACT 

We have mapped the subminimetre wavelength continuum emission from the Subaru 
Deep Field (SDF) at 450 and 850 /im with the Submillimetre Common- User Bolometer 
Array (SCUBA) detector on the James Clerk MaxweU Telescope (JCMT). The near- 
IR image of the SDF is one of the deepest near-IR images available and contains four 
'hyper extremely red objects' (HEROs). These data allow us to test the connection 
between 'extremely red objects' (EROs) found in IR surveys and the population of 
bright submillimetre sources found with SCUBA. We present a weak measurement of 
the average flux of the four JT'-band selected HEROs of 1.15 (±0.46) mJy, which fails 
to support the hypothesis that HEROs should be bright SCUBA sources. Our data 
are consistent with the HEROs being objects with SEDs like that of Arp220 out to 
z ~ 1.7, however, the extinction in the HEROs must be about 1 magnitude greater 
in the J-band than is the case for Arp220 and they would need to be 1.7 times as 
luminous as Arp220. On the other hand, an evolutionary model of elliptical galaxies 
at z ~ 2-3 in a dusty starburst phase is also in agreement with the submillimetre 
data, as was originally proposed for the HEROs. 

Key words: submillimetre - cosmology: observations - galaxies: high-redshift - 
galaxies: evolution - galaxies: formation - galaxies: starburst 



1 INTRODUCTION 

Extremely Red Objects (E ROs) emerge at faint near- 
IR magnitu des (K > 18) jElston. Riek e fc Rieke ) I l988t 
[McCartlw, P ersson fc Westlll992HHu fc Ri dgwav 19941) and 
have become important in studying the hist ory of galaxy 
formation and star formatio n in the Universe dYahata et alJ 
I2OOOI: iDickinson et al.ll200(]ll . The ERO class contains ob- 
jects with very red optical to near-IR colours, typically 
I — K > A or R ~ K > 5 — 6. This extreme observed 
colour can be produced by old stellar populations at 
2: ~ 1 or young dust-enshrouded starburst galaxies at 
high redshifts. Spectroscopic observations of EROs have 
revealed that about half of the population is comprised 
of very dusty galaxies undergoing their initial burst of 
star formation, while the remaining EROs are evolved stel- 
lar populations (|DunloB_et_alJ [1 9961: iGraham fc Devlll996l: 
ICohen et alJll999l : ICimatti et alJll999l) ! 

Strong submillimetre emission from an ERO is a n in- 
dication of starburst activity (see IWebb et al.l i2004ft '). In 
fact, many of the most luminous submillimetre sources dis- 
covered (with typical fluxes of 5 — 20 mJy at 850 ^m) have 
been identified as counterparts to EROs (e.g. CUDSS14A 



JOear et alJl200ol). SMM J00266-H7 08 iFraver et al."2000h. 
W-MMDll JChapman et all 120001) SMM J09429-H4658, 
and SMM J04431-H0210 JSmafl et al]|l999^ ). Because of the 
high frequency with which very bright submill imetre sources 
are identified with EROs, ISmail et al.l il999() suggest that 
SCUBA sources comprise the majority of the reddest EROs. 
This would imply that a Tf-band source which is very faint 
or undetected in optical data (i.e. J or /-band) would likely 
be a bright submillimetre source. Here, we test if the most 
extremely red objects found to date are strong submillimetre 
emitters. 

The Subaru Telescope has been used to produce deep 
images at J (A = 1.16-1.32 /xm) and K' (A = 1.96-2.30 /xm) 
of a 'blank' 2 arcmin x 2 arcmin field near the north 
Galactic pole with 5 a dete ction limits of 25. 1 and 23.5 in 
J and K' , respectively iMaihara et ablboOll) . Among the 
roughly 350 galaxies detected in the K'-hand in the Sub- 
aru Deep Field (SDF), there are 4 objects with extreme 
J — K' colours iMaihara et alj200J) . These 4 objects have a 
1 /error^ weighted average K' magnitude and J — K' colour 
of 21.6 and 3.0, respectively (see Table0. 

Among the EROs which are faintest in the K'- 
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ID Position (2000.0) K' J - K' 

SDFl 13''24™22!38 +27°29'49y5 20.91 (± 0.05) 2.97 {± 0.14) 

SDF2 13''24™22!39 +27°29'0iy9 22.03 (± 0.09) 3.65 {± 0.40) 

SDF3 13'^24™2m6 +27°29'0iy9 21.99 (± 0.05) 2.81 (± 0.20) 

SDF4 13'^24™22';84 +27°30'08y4 22.31 (± 0.14) 4.12 (± 1.04) 

Table 1. The positions, K' magnitudes, and colours of tlie HEROs from lMaihara et aljl200ll) . Astromctry is derived from the coordinates 
of an HST guide star in the flanking field and the estimated positional accuracy is ±0'.'15. 



band {K' > 22), an increasing fraction are extraordinar- 
ily red , even comp a red t o other EROs, with J — K' > 
3 or 4. iTotani et al] J200ll) call these objects Hyper-EROs 
(HEROs) and, in detailed modelling, find that they are too 
red to be passively evolving elliptical galaxies. The HEROs 
are most likely to be either very dusty galaxies which formed 
at redshifts of z ~ 4 — 7 and are still undergoing rapid star 
formation when seen near a 3, or they are clean 'Lyman- 
dropout' types of galaxies seen at extraordinary redshifts of 
a « 8 — 10. 2 mjy SCUBA sour ces have a surface density of 
3 X lO^deg"^ (iBorvs et al.l2003^ . which is similar to the SDF 
HERO surface density of « 3 x lO^deg"^. If HEROs com- 
prise the bulk of the SCUBA population, then their number 
density implies a brightness of 2 mJy. A measurement of the 
submillimetre flux of these HEROs would allow us to distin- 
guish between these models. 



One of the bright SDF HEROs may be an appar- 
ent interacting pai r of galaxies (see fig. 13 (SDFl) of 
iMaihara et aljl200ll where both components are clearly un- 
detected in the J-band), which could be exciting, since 
mergers have been known to play a role in models of ul- 
traluminous infrared galaxies (ULIRGs), which are proba- 
bly related to SCUBA sources. Us ing Poisson statistics and 
N ~ lO'^deg"^ down to 23 jMaihara et alJl200ll) . we 
estimate that there is only a 5 per cent probability that two 
galaxies would appear 1.5 arcseconds apart by chance (ignor- 
ing clustering). And with the previous result, we determine 
that there is an 80 per cent chance of not seeing one galaxy 
pair (within 1.5 arcsec) among four galaxies. However, for 
a HERO surface density of 4.5 x 10'' deg~^, the chance that 
two HEROs in the SDF would lie within 1.5 arcsec of each 
other at random is 0.25 per cent. It is therefore quite un- 
likely to see a pair of HEROs 1.5 arcsec apart if they were 
not interacting. One pair of HEROs in the SDF corresponds 
to 900 deg^'^. How bright do we expect a pair of HEROs to 
be in the submillimetre? If they represent a large fraction 
of bright submillimetre sources, then their surface density 
suggests a flux d ensity of around 10 — 20 mJy (see for e.g. 
iBorvs et allbOOSl) . 



The paper will be organized as follows. In Section |5| 
we discuss the observations, data reduction and fluxes. The 
SED galaxy templates and the expected 850 /im flux ver- 
sus redshift relation are discussed in Section ^ Section 2] 
presents the results and discussion, and Section |K| states the 
conclusions. Throughout this paper we assume standard cos- 
mological parameter values of Qa = 0.7, = 0.3 and 
Ho = 75kms-^Mpc-^ 



2 OBSERVATIONS, DATA REDUCTION AND 
FLUXES 

The SDF was observed with a resolution of 14.8 arcsec and 
7.5 arcsec, a t 850 and 450 Atm, res pectively, with the SCUBA 
instrument jHoUand et alJll999l) on the 15-m JCMT atop 
Mauna Kea in Hawaii. Observations were carried out dur- 
ing May of 2001, May 2002, February 2003, and March 
2003. The atmospheric zenith opacity at 225 GHz was mon- 
itored with the Caltech Submillimetre Observatory (CSO) 
tau monitor in 2001 and 2002 and with the JCMT water 
vapour monitor in 2003. 

One shift of 'jiggle mapping' was performed in May 2001 
in unstable medium-grade weather, where the 850 /im atmo- 
spheric zenith opacity ranged from 0.26 to 0.44. Jiggle map- 
ping, centred on an RA and Dec (J2000) of 13''24™21':20 
and -|-27°29'25'.'0, was performed using offset pointing cen- 
tres and the secondary mirror was chopped at a standard 
frequency of 7.8125 Hz in RA/Dec coordinates, with chop 
throws of 30 and 40 arcsec. Additional jiggle map data were 
obtained in March 2003 in dry weather, using chops in az- 
imuth in order to reduce the effect of rapid sky variations 
with alternating position angles of and 90 and chop thows 
of 30 and 40 arcsec. Pointing checks were performed hourly 
on blazars and planets and offsets from the pointing centre 
were typically 2-3 arcseconds. 

In addition, 9-point photometry data were obtained in 
May 2002, February 2003, and March 2003, using the cen- 
tral bolometer on the SCUBA array for each of the four 
HERO sources. Chop throws were selected to avoid the other 
sources in the fleld. In total, 13.5 and 16.5 hours of useful in- 
tegration time were obtained for the SDF at 450 and 850 /im 
respectively. 

We used SURF (SC UBA User Reduction Facility; 
IJenness fc Lightfootlll99 8^ together with locally developed 
code, written bv ]BOTv^ (2002). for data reduction. The re- 
duction steps were standard, except that the flux from the 
negative off-beams were then 'folded in' to increase the over- 
all sensitivity of the 850^m map by a factor of roughly ^3/2. 
We did not do this for the 450/im map, since the beam is not 
sufficiently Gaussian-shaped for this procedure to be simple. 

Calibration data were reduced in the same way as 
the SDF data. Due to the poor weather and time con- 
straints, few calibrators were observed during the photome- 
try run, and a standard ffux conversion factor was adopted 
for those data. Neglecting undersampled and noisier map 
edges, the maps have mean values of —0.05 (±0.03) mJy and 
0.3 (±0.2) mJy at 850 and 450 /^m, respectively. This is con- 
sistent with the expected map average of zero for differen- 
tial measurements. The rms values over the maps are about 
2.3 mJy and 11.4 mJy at 850 and 450 ^tm, respectively. 

The flux measurements for the HEROs also include 
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ID 


5450 




5850 95% limits 


1 




(mjy) 


(mJy) 


(mJy) 




SDFl 


1.5 (± 5.9) 


2.1 (± 0.9) 


0.4 < 5850 < 3.8 


u 


SDF2 


-7.8 (± 6.4) 


1.0 (± 0.9) 


0.0 < 5850 < 2.5 


V 1 


SDF3 


6.1 (± 6.5) 


-0.22 (± 0.9) 


■5850 < 1-6 


o 


SDF4 


1.7 (± 5.8) 


1.8 (± 1.0) 


0.1 < 5850 < 3.5 


cri 

(D 


Mean 


0.42 (± 3.06) 


1.15 (±0.46) 


0.3 < 5850 < 2.0 


Li_ 



Table 2. The measured 450 and 850 /xm flux densities and errors 
for the four SDF HERO sources. The fourth column lists the 95 
per cent Bayesian lower and upper limits to the 850 /im flux. The 
bottom row lists the error-weighted mean of the 850 /.tm fluxes 
and the 95 per cent Bayesian lower and upper confidence limit to 
this mean flux. The 95 per cent Bayesian upper confidence limit 
to the 450 /.tm flux is < 6.3 mJy. 



the photometry observations, which are treated as under- 
sampled maps and folded in with the jiggle map data. 
Since our targets are unresolved by the JCMT, we sim- 
ply measure the flux on the beam-convolved SCUBA map 
at the Tf'-selected position of each HERO. We then use 
the noise map to find the noise associated with each 
pixel, ignoring the calibration uncertainty, which is unim- 
portant at these low signal-to-noise ratios. The galac- 
tic cirrus contribution is 0.9MJys r~^ (extracted from the 
ISchlegel. Finkbeiner fc Davij|l993 map) which is much less 
than the confusion limit of the map. 

The results are summarized in Table |5| We fail to sig- 
nificantly detect any of the HEROs and therefore conclude 
that they are not luminous SCUBA sources. Since we can- 
not claim detection of any of our faint sources, Bayesian 
95 per cent upper confidence limit flux estimates are calcu- 
lated for each source by integrating over the non-negative 
flux regions of a Gaussian probability function. Combin- 
ing the results from the four sources allows us to obtain 
an average object flux for our sample of HEROs. The error- 
weighted average flux density for the whole sample at 850 /im 
is 1.2 (±0.5) mJy, and has a 95 per cent upper confidence 
limit of 2.0 mJy. The average fiux density at 450 /im is 
0.4(±3.1)mJy. 

If HEROs comprise a large fraction of the submillime- 
tre source counts of Borys e t al. (2003), the measured mean 
brightness of 1.15 (±0.46) mJy would correspond to sources 
with a mean density of 6.7 (±3) sources per 2arcmin x 
2 arcmin area of sky, consistent with the 4 sources found 
in the Subaru Deep Field. 



3 THE MODEL 

In order to interpret our results, we need a realistic model 
to describe expectations for the submillimetre properties of 
the HEROs. We use a phenomenological model, based on 
scaling the properties of known prototypes. We utilize the 
850 /^m fiuxes, K' magnitudes and near-IR colours, to fit 
a three parameter model to a redshifted SED model. The 
model parameters include the overall luminosity, the amount 
of dust reddening in the IR, and the redshift. 



<D 1 




10' 10^ 10-- 
Observed Wavelength [/xm] 

Figure 1. The model spectral energy dist ribution of a nearby 
ULIRG: the iLagache. Dole fc Pugetj ^2002^ starburst template 
(dashed line represents number 30 of 41 possible templates) fit 
Skt z = 0.018 to Arp220 photometric data (solid line is an ex- 
tinction corrected version), which are marked with the plus signs. 
The slope of the drop at < 1 /jm has been extrapolated to shorter 
wavelengths. 



3.1 ULIRG model 

We created a template SED using public photometric data 
available for Arp 220 from the NASA/IPAC Extragalactic 
Database (NED). Arp 220 is the most luminous object in 
the nearby Universe {z — 0.018) and is a well-studied exam- 
ple of a ULIRG. Higher redshift sources have more sparsely 
sampled SEDs, and so in using Arp220, we start at least on 
firm ground. 

We then use the iLaeache. Dole fc Puge3 ^2002^ UV- 
to-radio SED models of a typical ULIRG to interpo- 
late the Arp220 photometric points (see Fig. 0. These 
templates evolve in luminosity (i.e. become hotter) and 
cover a wide range of luminosity levels from L = 10^ — 
IO^^Lq. We employ their model number 30. We increase 
the template's IR dust content (see details below) since the 
Lagachc. Dole fc P ugct (2002) starburst template is not red 
enough to fit the Arp220 photometric points. In addition, in- 
corporating a dust model provides a mechanism for varying 
the amount of dust extinction to fit the HERO colours. 

We redden the original template by multiplying it by 
an extinction function which varies with wavelength: 



/= 1 



1.23 



(1) 



where the free parameter, x, is adjusted to a value of 
X = 2.5 in order to match the NED data points, and 
A is wavelength in /im. The above extinction law was 
constructed using a power-law fit to the relative ex- 
tinction in the Landolt V, R, I and UKIRT J, H, 
K and L' bandpasses versus effective wavelength from 
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1 0-2° I 

10-1 10" 10' 102 10^ lO-* 10^ 

Observed Wavelength [/im] 

Figure 2. The model sp ectral energy distribution of a nearby 
normal spiral galaxy: the iLagache. Dole fc Puge3 ^2002^ normal 
galaxy template fit at ^ = 0.002699 to NGC3938 photometric 
data, which are marked with the plus signs. Again, the model 
SED has been extrapolated to shorter wavelengths. 

ISchlegel. Finkbeiner fc Davi^ (|l99?j ). They employ a 'dif- 
fuse IS M' mean value of Ry = 3.1 fo r the extinction 
laws o f ICardelli. Clayton fc Mathi3 (Il989l) and lO'Donnelll 
1^9^. We note that the difference between the extinc- 
tion curves of the Milky Way, the Magellanic Clouds and 
starburst galaxies is almost negligible at wavelengths longer 
than ~ 0.26 Mm (see [ C alzetti^ Kinney fc Storchi-BergmannI 
Il994l: ICardelli. Clavton fc Mathij Il989^ ■ Hence variations 
are unimportant for z < ?> galaxies (where the J-band cor- 
responds to rest wavelengths of < 0.31 ^m). 

3.2 Spiral Galaxy Model 

For comparison, we also created a normal spiral galaxy 
SED using NED pubUc photometric data for NGC3938. 
NGC3938 is a well-studied, multiple-armed, early luminos- 
ity class Sc(s) face-on g a laxy 1 0.8 Mpc away. We utilize the 
iLaeache. Dole fc Puged (|2002) template of a normal cold 
spiral galaxy scaled evenly across all wavelengths with no 
addition of dust required (see Fig. |5J . 

3.3 Using the Templates 

We have constructed three-parameter models based on these 
SEDs, constrained only by the K' magnitude and the J — K' 
colour. We choose a redshift, redden the template to match 
the observed J — K' colour by varying x in Eqn. Q then 
boost or diminish the template achromatically to match the 
K' magnitude. 

We note that performing numerical integration over the 
Arp220 SED of Fig. ^reveals that the total power absorbed 
in reddening the starburst galaxy template is only about 3 — 



5 r < < < < \ < < < < \ < < < < \ < < < < : 
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0.0 0.5 1.0 1.5 2.0 

Redshift z 

Figure 3. The amount of dust extinction that we must add to 
our modified Arp220 SED (bold lines) and to our normal spiral 
galaxy SED for increasing redshift at the J (solid lines) and K' 
(dotted lines) bands to achieve the mean observed J — K' colour. 

5 per cent of the total galactic emission. The small increase 
in the far-lR luminosity can be neglected as the reddening 
is altered. In contrast. Fig. |5| reveals the energetics of the 
normal cold spiral galaxy, and it is clear that there is just 
as much power in the near-IR peak as in the far-IR peak. 
Therefore, the power that would be 'lost' due to reddening 
the spiral galaxy must be added back into the far-lR region 
of the spectrum in order to conserve energy. We adjust the 
FIR luminosity accordingly in order to account for the power 
absorbed in reddening the galaxy. 

We scale the template brightness uniformly across each 
spectrum to fit the inferred K' flux, as determined from the 
observed K' magnitude using: 

= lO-'^'^'^Fo, (2) 

where m is the magnitude and Fq is the zero-point flux which 
is 718.9 Jy for the i^'-band, and was extracted from the 
interpolation of the broad-band flux of Vega. 

After fitting the template to the K' magnitude and 
J — K' colour, we then read off the observed 850 ^m flux. 
We repeat this process for a range of redshifts (see Fig. 01. 
Added dust in the J and K' bands versus redshift is plotted 
in Fig. |3] We plot the scale factor versus redshift in Fig. 2] 



4 RESULTS AND DISCUSSION 

We estimate a redshift for each HERO using the recipe that 
follows. At each redshift, we redden the template to produce 
the observed J — K' colour of each HERO. Next, we scale the 
whole SED in order to match the observed K' magnitude. 
We then test whether the predicted 850 fim flux is consistent 
with the observed value. Figs.|3|3]and|S|can be read to find 
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Redshlft z 

Figure 4. A plot of the brightness in the i^T'-band (giving the cor- 
rect K' magnitude) compared to the brightness of Arp220 (lower 
bold line) and a normal spiral galaxy (upper line) after adding 
in the appropriate amount of extra dust to obtain the observed 
J — K' colour, plotted versus redshift for our moan HERO. 
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Figure 5. Expected 850 /im flux density versus redshift for our 
mean HERO assuming a modified Arp220 SED (upper bold line) 
and a normal spiral galaxy (lower line) . Asterisks represent where 
we read off a redshift, based on the mean and upper limit fluxes. 



ID Upper Limits 





2ULIRG 


■^spiral 


SDFl 


< 1.7 


< 1.8 


SDF2 


< 1.6 


< 1.7 


SDF3 


< 1.8 


< 2.0 


SDF4 


< 1.6 


< 1.8 


Mean 


< 1.7 


< 1.8 



Table 3. Model-estimated redshifts (for the ULIRG and the nor- 
mal spiral galaxy) based on the 850 ^m upper limit flux of each 
source. The bottom row lists the model-estimated redshift based 
on the upper limit flux of our mean HERO. 

the properties of a ULIRG or a spiral galaxy corresponding 
to our measurements. 



4.1 Could the HEROs be Starburst Galaxies? 

Using Fig. |3 we see that an object with an SED like that 
of Arp220 and properties corresponding to the mean of 
our measurements in Table |21 must be located a.t z = 
1.44 (zbO.lO). From Fig. 2] one sees that such a reddened 
object at z — 1.44 must be 0.9 (±0.2) times as luminous as 
Arp220 to match the observed i^^'-band flux, and from Fig.|21 
the amount of reddening in K'^bseived is 0.4 magnitudes com- 
pared to Arp220. Our observed 95 per cent upper confidence 
limit for the four HEROs of Sgso < 2.0 mjy leads in the same 
way to z < 1.7 (see Table EJ, L^{HERO)/-t'..{ULiRG) < 1-7, 
and AK > 0.4 magnitudes of extinction. This is consistent 
with the median re dshift estimates f o r SCU BA sources from 
iDunlotj ll200ll) and lChapman et al.l J2003h . but it is at the 
low redshift end. If the HEROs have SEDs Uke Arp220 and 
arc located at the mean redshift reported bv lChapman et all 
(2003) ot z — 2.4 and possess the observed J and i<''-band 
fluxes, one would expect them to be 12mJy sources, which 
they are not. 

In principle the 95 per cent upper confidence limit of 
<S'45o < 6.3 mJy provides a redshift upper limit ot z < 1.44, 
which is tighter than the limit obtained from the 850 j-im 
data. However, the value of this limit relies on the precise 
understanding of the 450 variance and we therefore do 
not rely on it. 

We have tested our Arp220 model by comparing it to 
those very few objects for which IR colour, if-band mag- 
nitude, 850 /xm flux and redshift are all known. We use an 
ERO and a VRO (very red object) from the Hubble Deep 
Field (HDF) (Borys ot al. 2003) wit h measured spectro - 
scopic redshifts. We also use HDF850.1 jHuehes et al.ll998l) . 
a SCUBA- bright IR-faint ERO with a secure photomet- 
ric redshift llDunlop et al.ll2004ft. We al so include HRIO (or 
ERO J164502-H4626.4. iDev etaPligggl) . the first ERO dis- 
covered to be associated with a ULIRG, and FNl 40 from the 
Far InfraRed Bac kground (FIRBAGK) survey (lDolell200d. 
ISaiina et aDl2003l) . See Table H for the IR colours, submil- 
limetre fluxes, spectroscopic redshifts and model-predicted 
redshifts. The mean difference between the spectroscopic 
redshifts and the model-estimated redshifts (^modei — ^spoc) is 
-1-0.06 with a spread of ±a = 0.6. We are therefore confident 
that this simple model is suitable for estimating redshifts of 
the reddest population of objects when the IR colour and the 
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850 /xm flux are known. The model possesses an advantage 
over more complicated multi-parameter models in that it fits 
known sources well, despite its single-parameter simplicity. 

Since the HEROs are extreme examples of EROs, it is 
useful to consider how the resul ts would change by fi tting 
an ERO SED to the data. HRIO jHu fc Ridgwavll^94^ is an 
archetypal ULIRG-ERO association and is merely a distant 
clone of Arp220 ( z = 1.44) and 3.8 times more luminous 
iElbaz et al.l l2002^. We construct a template of HRIO and 
find nearly identical constraints on redshift an d lumin osity. 

We treat the mo del of iTotani et al] ll200lll (see 
iTotani fc Takeuchill200 j for details of the model) of proto- 
elliptical galaxies which formed at 2 = 3 and are observed 
at 2 = 2.3 as a template for emission. We apply Eqn. 
to redden the template and we find that achieving a colour 

01 J ~ K' = 4.12 requires 1.80 magnitudes of extinction in 
J and 0.57 magnitudes of extinction in K' . The resulting 
850 /iui fi ux is 0.5 mjy . This is close to the 1.0 mjy predic- 
tion from ITotani et aO i200 jl and consistent with our mea- 
sured average flux of 1.15 (±0.46) mJy. However, there are 
^ 2 X 10"* 0.5 mJy SC UBA sources per d eg^, according to 
the number counts of iBorvs et alj (|20o3). Recall that the 
surface density of HEROs in the SDF is ~ 3 x 10^ deg~^ 
which corresponds to the surface density of 2 mJy SCUBA 
sources. If HEROs are typical elliptical galaxies then they 
only comprise a small fraction of SCUBA-detected sources. 

4.2 Could the HEROs be Normal Spiral Galaxies? 

We now turn to the normal spiral galaxy SED. When we 
force our template SED to match the 95 per cent upper 
limit to the 850 /xm flux of 2.0 mJy, constrained by a faint 
K' magnitude and very red J — K' colour, we estimate a 
redshift oi z < 1.8 (see Fig.|^. 

At 2: = 1.8, the SED requires a factor of ~ 32 in luminos- 
ity boosting in order to reproduce the correct K' magnitude 
for the HERO sample (see Fig. |1J. At 2 = 1.8, we require 
r; 0.4 magnitudes of extinction in K' and ~ 1.2 magnitudes 
of extinction in J (see Fig.|3J. 

Because of the large luminosity boosting required, 
NGC3938 is not a reasonable analogue for HEROs if they are 
near to the 850 fim flux upper limit. However, if the HEROs 
are normal galaxies at a redshift near to 0.6 and reddened by 

2 magnitudes in J-band, they would produce the observed 
J and /f'-band fluxes with no luminosity boosting required. 
The 850 ^im flux would be 0.06 mJy. This is below the 95 
per cent lower confidence limit of Sg5o > 0.3 mJy, although 
the value of this constraint depends on the reliability of our 
estimate of the variance in the data. 



4.3 Otiier Possibilities 

Can the HEROs be evolved elliptical galaxies at higher red- 
shifts? We know that the reddest and faintest EROs at 
known redshifts are not elliptical galaxies. Therefore, the 
HEROs are most likely not well-described by high-redshift 
passively evolving elliptical galaxies that are red due to old 
stellar populations. 

If the HEROs are J-band d ropout galaxies, they must 
lie at redshifts of 2 8 - 10. ITotani et all ll200ll) state 
that the UV luminosity inferred from the K' magnitudes of 



the HEROs implies star formation rates of > 140 M© yr~ , 
which sets a total mass of M > 10^ Mq for each system. 
They find that the calculated surface density of this popula- 
tion of massive objects at such early times (in the standard 
theory of CDM hierarchical structure formation) is much 
too small compared with the observed number density of 
HEROs in the SDF. Thus, there is no reason to suspect that 
these objects are J-band dropout galaxies at extraordinary 
redshifts. 



5 CONCLUSIONS 

We observed four objects with extreme J — K' colours 
(HEROs) with SCUBA on the JCMT. We obtained no sta- 
tistically significant detections of the individual objects, but 
we obtained a marginally significant detection of the objects 
as a group, finding a mean 850 /im fiux of 1.15 (±0.46) mJy 
or a combined 95 per cent upper confidence limit of 2.0 mJy. 
This is no more than a hint that they are dusty galaxies 
emitting a small amount of submillimetre fiux. 

HEROs could have SEDs hke that of Arp220, but the 
ones we examined must reside at the small redshift end of 
submillimetre redshift distribution, with an average redshift 
of 2 « 1.4. These objects would have rest frame luminosities 
comparable to that of Arp220 and would be only slightly 
reddened. 

HEROs could alternately be the progenitors of ellip- 
tical galaxies at t he end of a dusty starburst phase, as 
ITotani et alJ i200ll) suggest. In this case, they would be at 
the dim end of the submillimetre fiux range that our mea- 
surements allow, and they do not comprise the majority of 
850 /im sources at that flux level. 

HEROs are not likely to be ordinary elliptical or spiral 
galaxies at any redshift. 

In order to conflrm that HEROs are bona flde SCUBA 
sources (which requires detecting individual objects with 
SCUBA), a much larger survey would need to be conducted. 
Assuming that the HEROs comprise a subset of the sub- 
millimetre source counts, a near-IR survey of 15 arcmin^ 
to the depth of the Subaru Deep Field survey should re- 
veal a HERO corresponding to an 850 fim source brighter 
than 5mJy with 95 per cent confidence. Our model indi- 
cates that sources brighter than 5 mJy would be at redshifts 
above 2 = 2.0, which is consistent with the distribution of 
bright submillimetre sources for which redshifts are known. 

SCUBA galaxies are very often identified with EROs 
and we have tested if the inverse statement also holds true. 
Our results clearly demonstrate that measuring the submil- 
limetre fiux of HEROs will not necessarily select very bright 
SCUBA sources. SCUBA galaxies may just comprise a sub- 
set of a more diverse population of EROs. 
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ID 


Optical 


^850 


^spcc 


^mo del 




/ 


I-K 


(mJy) 






SMMJ123628+621048 


22.5 


4.0 


4.4 ± 1.2 


1.013 


1.2 ± 0.3 


SMMJ123635+621239 


22.3 


3.5 


3.0 ±0.8 


1.219 


1.3 ± 0.3 


SMMJ123652+621227 


28.7 


5.2 


7.0 ±0.5 


4.1 ± 0.5 


3.4 ±0.1 


SMMJ164502+4626.4 


24.9 


6.5 


4.9 ±0.7 


1.44 


0.82 ±0.1 


FIRBACK FNl 40 


23.96 


4.23 


6.3 ± 1.4 


0.449 


1.8 ± 0.2 



Table 4. Submillimetre sources with red IR colours and known redshifts. Model-predicted redshifts and errors (estimated using the 
errors in the 850 /im fluxes) are given in the last column. 
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